New isospin effects in central heavy-ion collisions at Fermi energies 
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Isospin effects on multifragmentation properties were studied thanks to nuclear collisions between 
different isotopes of xenon beams and tin targets. It is shown that, in central collisions leading 
to multifragmentation, the mean number of fragments and their mean kinetic energy increase with 
the neutron-richness of the total system. Comparisons with a stochastic transport model allow to 
attribute the multiplicity increase to the multifragmentation stage, before secondary decay. The 
total charge bound in fragments is proposed as an alternate variable to quantify preequilibrium 
emission and to investigate symmetry energy effects. 
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I. INTRODUCTION 

One of the present motivations for investigating heavy- 
ion collisions at intermediate energy is the improvement 
of the knowledge of the Equation of State (EOS) for nu- 
clear matter. More specifically the formulation of an 
adequate symmetry term is required to progress, and 
the density dependence of this term both at sub- and 
supra-normal density is still debated. Further experi- 
mental constraints are clearly necessary. Within the next 
decade physicists expect the advent of new heavy-ion ac- 
celerators, providing high-intensity exotic beams, in or- 
der to study reactions covering a broad range of isospin 
(N/Z) ratios; in the meantime information on isospin 
effects can be obtained thanks to a judicious choice of 
projectile-target couples. Comparisons of several cal- 
culated and experimentally measured isospin-dependent 
variables already provided some hints on the symme- 



* Current address: PTB Brauncshweig, Bundesallee 100, 38116 
Braunschweig, Germany 

t rivet@ipno.in2p3.fr 

* Present address: S2HEP (EA4148), UCBL/ENSL, Universite dc 
Lyon, Villeurbanne, France 

§ http://indra.in2p3.fr/spip/ 



try term, but the results are still highly model- and 
experiment-dependent [1-7] . 

In this line the INDRA collaboration has studied colli- 
sions between 124,136 Xe projectiles, with an incoming en- 
ergy of 32 and 45 A MeV, and 112 < 124 Sn targets. We show 
in this paper that, in central collisions, the mean charged- 
product and fragment (Z>5) multiplicities, and mean 
fragment kinetic energies, depend on the isospin of the 
total system. The commonly accepted reaction scenario 
for central heavy-ion collisions around Fermi energy, vali- 
dated by stochastic transport models, is the following [8]: 
there is firstly a compression phase, the strength of which 
depends on the masses and mass asymmetry of the in- 
cident partners. This stage is followed by an expansion 
phase accompanied by the emission of fast preequilibrium 
particles. The diluted system enters the spinodal region 
of the phase diagram and eventually breaks into several 
excited fragments, light particles and nucleons (multi- 
fragmentation); when the configuration is frozen, namely 
when the nuclear interaction between fragments becomes 
negligible, one speaks of the "freeze-out" stage. Finally 
the charged products move apart, further accelerated by 
the Coulomb force, while losing their excitation energy 
through evaporation. 

In this framework we shall test, through comparisons 
with a stochastic transport model, the origin of the 



isospin dependence observed for the aforementioned vari- 
ables. 



II. EXPERIMENT 

Beams of 32 and 45 A MeV 124 < 136 Xe accelerated by the 
Grand Accelerateur National d'lons Lourds (GANIL) im- 
pinged on 530 /jg.cm~ 2 112 ' 124 Sn targets. Beam intensity 
was about 3-5 x 10 7 particles per second to avoid event 
pile-up. Charged products were detected and identified 
with the Ait INDRA array [9] , which comprises 336 two- 
or three-member telescopes arranged in a cylindrical ge- 
ometry around the beam axis. The array was upgraded 
as follows [10]: 12 (300 ^m)silicon-(14 cm)CsI(Tl) tele- 
scopes replaced the phoswiches of the first ring, and 7 of 
the 300 fim silicon wafers were replaced by 150 /im wafers, 
with an increased amplifier gain, in order to identify iso- 
topes up to nitrogen. In the standard silicon-CsI(Tl) 
telescopes (2-45°) elements were identified within one 
charge unit; in addition between 14° and 45° these tele- 
scopes provided isotopic identification up to carbon. The 
same Z-resolution was obtained in the first stage ioniza- 
tion chamber-silicon or ionization chamber-CsI(Tl)) tele- 
scopes up to Z=20. For low energy heavier nuclei, Z are 
known with an uncertainty of 1-2 charge units [11]. In 
this paper, the term "fragments" refers to charged prod- 
ucts with a charge Z ^ 5 whereas light charged particles 
(lep) stands for particles with a charge Z ^ 2. 

Most of data taking was done with an on-line trigger 
requiring at least four fired INDRA telescopes. For a 
better overview of the reactions, some inclusive measure- 
ments (one telescope fired) were also performed. 

Absolute cross sections were derived from the mea- 
sured target thicknesses, the counting of ions collected in 
the Faraday cup located at the end of the beam line and 
the acquisition dead time. The charge of ions reaching 
the cup was obtained using the formulae of reference [12]. 
The error on cross sections is estimated to 15%, while rel- 
ative errors between two systems are lower, around 3%. 

Table I summarizes the studied systems and the isospin 
contents of the total systems. Note that two of the sys- 
tems have the same N/Z, allowing to study entrance 
channel mass-asymmetry effects. We also indicate the 
measured inclusive cross sections (trigger multiplicity 
>1). In the off-line analysis we required at least one 
identified charged product and rejected events formed of 
a single charged product with an atomic number close to 
the projectile one, in order to exclude elastic scattering. 
The cross sections so obtained are close to the reaction 
cross sections calculated with the systematics of Kox [13] 
at 32 AMeV whereas they are smaller by about 1 barn 
at 45 AMeV. 

Figure 1 shows the charged-product multiplicity dis- 
tributions of identified nuclei at 32 A MeV, for data with 
an on-line trigger multiplicity M >1. For low multi- 
plicities the cross sections are higher for the neutron- 
rich system; the trend is reversed for multiplicities larger 



Projectile Target (N/Z) E p 



syst. (A MeV) (b) (A MeV) (b) 



124 



Xe 



112 Sn 1.269 

124 Xe 124 Sn 1.385 

136 Xe 112 Sn 13g5 

136 Ve 124 Sn 1.50 



32 5.16 45 4.73 

32 5.14 

32 5.90 45 

32 5.81 45 5.02 



TABLE I. Available systems and measured cross sections. 
Systematic errors on cross sections are estimated to 15% 
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FIG. 1. (color on line) Charged-product multiplicity distri- 
butions for the different systems at 32 A MeV. Error bars are 
statistical. 



than 20, and the distribution extends to higher values for 
the neutron-poor system. The two systems 124 Xe+ 124 Sn 
and 136 Xe+ 112 Sn display intermediate values. A first 
isospin effect appears already on raw variables: higher 
charged products multiplicities are observed when the 
system contains less neutrons; it is indeed expected that 
a neutron-rich system will preferentially emit neutrons to 
the detriment of light charged particles, which dominate 
in the charged product multiplicity. 



III. COMPACT SHAPE EVENTS FROM 
CENTRAL COLLISIONS 

In the following we shall report on data obtained with 
an on-line trigger M >4 for which we got a very high 
statistics. As in previous studies [14-17] we select quasi- 
complete events by requiring that the sum of the charges 
of the detected products, Z to t, be at least equal to 80. 
Due to the response of INDRA to the kinematics of 
quasi-symmetric collisions such as Xe+Sn, this choice 
favourizes de facto central collisions. We then isolate 
compact shape events (quasi-fusion) through the addi- 
tional condition that the flow angle (Ofiow) be greater 
than 60°. Let us recall that Qn w characterizes the 
main direction of matter emission in the reaction center 
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FIG. 2. (color on line) Evolution of the light charged particle 
and fragment multiplicities versus the N/Z of the total sys- 
tems, at the two energies. Error bars are smaller than the 
size of the symbols. 



of mass and is determined by the kinetic energy flow- 
tensor calculated from fragment momenta. More than 
10 5 events were selected in all cases, corresponding to 
measured cross-sections of 30-40 mb at 32 A MeV, and 
20-25 mb at 45 AMeV. These values are slightly larger 
than those previously measured for the 129 Xe+ nat Sn sys- 
tem [16]. The total cross-section for quasi- fusion, tak- 
ing into account detection efficiency and selection biases, 
is estimated to reach ^200-250 mb at 32 A MeV, and 
^170 mb at 45 A MeV (the difference between estimated 
and measured cross-sections is due to the limits on 9fi ow 
and on Z to t imposed in the selection, which cut part of 
events that can be attributed to the quasi-fusion mecha- 
nism [18]). 



A. Multiplicities 

Table II displays the different average multiplicities 
measured for the selected class of events, for all studied 
systems: M tot is the total charged product multiplicity, 
Mi cp and Mf rag refer to the light charged particle and 



fragment multiplicities, which are shown in fig. 2. Results 
from the reaction 129 Xe+ na *Sn have been added [16]: 
the composite system is the same as those formed by 
136 Xe+ 112 Sn and 124 Xe+ 124 Sn (the average mass num- 
ber of raa *Sn is 119). As in the unsorted data, at each 
beam energy, the average multiplicities of charged prod- 
ucts and of light charged particles decrease for increasing 
N/Z of the total system. Conversely the fragment multi- 
plicity increases. The three systems with the same N/Z 
show very close mean lep multiplicity values, which in- 
dicates that the entrance channel mass asymmetry has 
a small influence. The same is not true for mean frag- 
ment multiplicities, whereas their standard deviations are 
equal. While we may suspect some difference in the IN- 
DRA response between the present data and those of [16] , 
such an explanation does not hold for 124 Xe+ 124 Sn and 
136 Xe+ 112 Sn at 32 AMeV. In that case we might en- 
visage some physical reason, which requires further in- 
vestigation. Previous studies [19, 20] already indicated 
that more fragments are emitted for neutron-rich sys- 
tems. The available energy per nucleon, at the same in- 
cident energy, is the same for all systems, within 1%, so 
the increased number of fragments cannot be explained 
by a larger available energy. Moreover Mf rag modestly 
increases between 32 and 45 A MeV: it was indeed shown 
in [16] that for the system 129 Xe+ na 'Sn Mf rag presents 
a maximum for an incident energy around 40 A MeV. 

In previous works the difference in fragment multiplic- 
ities was attributed to phase-space effects [20], because 
statistical calculations (EES [21]) reproduce the observa- 
tion, or to sequential decay effects [19]. In the present 
data the fragment multiplicity increases by 10%, both at 
32 and 45 AMeV, between the lightest and the heavi- 
est systems. This corresponds to the mass increase be- 
tween these systems, which may recall the scaling law ob- 
served in [22], expected if multifragmentation originates 
from volume instabilities (spinodal decomposition). Us- 
ing stochastic transport codes we can test whether the 
increased fragment multiplicity arises from the dynami- 
cal or from the secondary decay stage of the reaction. 



B. Fragment kinetic energies 

The average fragment kinetic energies (in the center 
of mass) provide another piece of information on possi- 
ble isospin effects. Figure 3 displays this observable for 
the different systems at the two incident energies, ex- 
cluding the largest fragment of each partition. The gen- 
eral aspect is the same as that previously observed for 
129 Xe+™ a *Sn [15, 17, 23], namely an almost linear rise of 
the average energy when the atomic number increases. 
Conversely the mean energy of the largest fragment (see 
one example in fig. 3) increases for small Z and then 
decreases. We confirm that, for a given Z, the mean 
energy is smaller when it corresponds to the largest frag- 
ment [17, 23]. The average kinetic energy for a given 
element is slightly larger (by 2-3 MeV) at the higher in- 
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TABLE II. Average measured multiplicities of charged products, M to t, of light charged particles (H, He), Mj CP) and of fragments, 
Mfmg for the different systems studied. Number in parentheses are the standard deviations of the corresponding distributions. 
Statistical errors on mean values are smaller than 0.01 for Mtot and Mi cp and Mj rag . Results on line 3 are from [16]. 



E/A=32 MeV 


E/A=45 MeV 


System 
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{Mfrag) 
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{M lcp ) 


{Mfrag) 


124 Xe+ 112 Sn 


25.12 (2.90) 


19.66 (3.24) 


4.11 (1.16) 
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FIG. 3. (color on line) Average fragment kinetic energy, in the reaction center of mass, vs the fragment atomic number for the 
selected events at 32 (left) and 45 (right) A MeV. The largest fragment of each partition has been removed. The average kinetic 
energy of the largest fragment is displayed for the 124 Xe+ 112 Sn at 32 A MeV. Open symbols for 136 Xe projectile, closed symbols 
for 124 Xe projectiles. Circles, squares and triangles: experimental data; stars: SMF simulations (cold fragments, filtered). 



cident energy. 

The present data provide new information: at a given 
incident energy the mean fragment kinetic energy de- 
pends on the isospin of the system; the heavier the mass 
of the system, the larger the fragment kinetic energy is. 
We have verified that the non-measurement of the masses 
of the heavy fragments is not responsible for the observed 
effect. Indeed the mean measured kinetic energies, in the 
laboratory, are larger for the neutron-rich system. The 
mean measured masses for Z=5-8 are also larger. 

The fragment kinetic energy is the superimposition of 
a disordered thermal motion and an ordered term. The 
former term comes from the temperature of the fragment- 
ing system at freeze-out; the latter one, radially directed, 
originates from the Coulomb force, which depends on the 
system and fragment charges, and an expansion term pro- 
portional to the fragment mass [24-26] . At a given avail- 
able energy the thermal term is expected to be identical 
for all systems. The same holds for the Coulomb term 
as the system charges remain very close even after pree- 



quilibrium emission (see next section). Moreover the in- 
crease of the fragment multiplicity with the system mass 
should slightly decrease the Coulomb part [27], which 
is not the observed effect. In this context, the increase 
of the fragment kinetic energy signals the larger masses 
of the primary fragments when the system initially con- 
tains more neutrons. This is true if the expansion energy 
per nucleon does not depend on the initial isospin of the 
system, at a given incident energy. Isospin differences be- 
tween systems thus seem to survive preequilibrium emis- 
sion [28]. 

IV. COMPARISON OF DATA WITH A 
STOCHASTIC MEAN FIELD MODEL 

We used the stochastic mean field (SMF) calculation 
described in ref. [29]. Simulations were performed for 
collisions corresponding to all systems (table I), up to 
a time equal to 300 fm/c. The isoscalar EOS is soft 
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(i£' oo =200 MeV), and two parameterisations of the po- 
tential part of the symmetry energy are used [30], an 
asystiff one linearly increasing with density while the 
asysoft form (SKM*) has a maximum around normal 
density. The collision term uses the free nucleon-nucleon 
cross-section, with its isospin, energy and angular de- 
pendences 1 . This version of the model is the same as the 
one used in ref. [30], in which information on the asy- 
stiffness of the EOS was derived from isospin diffusion in 
the Ni+Au system at 52 and 74 A MeV. There is a great 
interest in finding additional experimental constraints on 
the density dependence of the symmetry energy in the 
same theoretical framework. 

At 45 A MeV we observe the formation of a single 
source which subsequently breaks into several fragments. 
At 32 A MeV, conversely to the Brownian One Body 
(BOB) calculations [31] shown in [16, 17, 32], the systems 
do not multifragment in head-on collisions. However it 
should be noticed that while the BOB calculations were 
perfomed for a thermalizcd source mimicking the com- 
posite system at the moment of maximum compression, 
here the whole dynamical evolution of the system, from 
the beginning of the reaction, is simulated by the SMF 
approach. Hence the impact of pre-equilibrium effects on 
the dynamics could be different in the two treatments. 

Fragments are recognized by applying a coalescence 
procedure to the one-body density, connecting nearby 
cells in which the density is larger than a cut-off value, 
taken equal to po/5 ("liquid phase"). We have shown 
in [32] that, at 300 fm/c, the fragment multiplicity is 
independent of the exact value of the cut-off density. 
The remaining early emitted nucleons constitute the "gas 
phase". The fragment phase space configuration at 300 
fm/c is injected in the SIMON code [33] which performs 
the secondary decay during the propagation of all prod- 
ucts under the Coulomb field, thus preserving space-time 
correlations. Note that the fragment excitation energies 
at 300 fm/c (~ 3.3 A MeV) agree well with experimental 
determinations [15, 34]. 1000 events were run with SMF 
for each system, then each primary event was de-excited 
20 times. The cold products were finally filtered through 
a geometrical replica of the INDRA array described in the 
code Panforte [35]. Because in the present calculations 
free nucleons ( "gas phase" ) emitted along the dynamical 
evolution are not taken into account in the subsequent 
de-excitation and Coulomb propagation steps, the selec- 
tion on the total detected charge cannot be used. More- 
over only fragment properties can be compared to the 
experimental values. 
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FIG. 4. (color on line) Fragment multiplicity distribu- 
tions, normalised to the number of events, for 124 Xe+ 112 Sn 
at 45 AMeV: experimental (circles), primary fragments 
(squares), cold fragments before (stars) and after (triangles) 
filtering. The asystiff EOS is used. Statistical errors are 
smaller than the symbol sizes. 



A. Fragment multiplicities 

We firstly verified the pertinence of the simulation 
by comparing calculated and measured distributions for 
some observables. We found that the charge distribu- 
tions are reasonably reproduced [36]. We show in Fig. 4 
an example of fragment multiplicity distributions: for 
the 124 Xe+ 112 Sn at 45 ,4 MeV the measured experimen- 
tal distribution is plotted together with those of primary 
fragments, cold fragments, and the filtered distribution. 
The primary and cold fragment distributions are very 
similar, both the width and the average value of the 
latter being only slightly smaller than those of the for- 
mer. Let us recall that we consider fragments starting 
from boron (Z=5). Two reasons support this choice in 
the simulations: first the yield of primary fragments is 
less sensitive to the fragment formation process, which 
comes from spinodal decomposition in the SMF model 
and strongly disfavours the formation of very small frag- 
ments [37]. Second these fragments are essentially rem- 
nants of primary hot fragments, whereas final Li and Be 
isotopes are also populated by secondary decay 2 . The 
filtered distribution shows a good agreement with the 
experimental one. 

The full set of data is represented on Fig. 5. Average 
multiplicity values are reported as a function of the N/Z 
of the sources, taken at t=120 fm/c (see next subsection). 
The experimental points are located at the average N/Z 
of the asystiff and asysoft sources. We first observe that 
the multiplicity of primary fragments (squares) increases 
with the isospin of the multifragmenting source, which 



2 The multiplicity of fragments with Z>3 is, for all systems, larger 
1 with an upper limit of 50 mb, to repress spurious low-energy after de-excitation than at 300 fm/c, which is not the case for 

collisions that of fragments with Z>5 
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FIG. 5. (color on line) Average fragment multiplicities, for the 
different systems at 45 AMeV: experimental values (circles, 
the line is to guide the eye), calculated values for primary 
fragments (squares), cold fragments before (stars) and after 
(triangles) filtering, are plotted vs the source N/Z. Asystiff 
EOS: open symbols, asysoft EOS: filled symbols. In all cases 
statistical error bars are smaller than the symbol sizes. 

allows to state that the experimental trend is essentially 
due to the preequilibrium stage of the collision. It would 
support the scaling law of multiplicity with system mass 
expected from a spinodal decomposition. The cooling 
of the fragments only slightly decreases the multiplici- 
ties, and more when the projectile is 124 Xe; indeed frag- 
ments from the less neutron-rich system are expected to 
evaporate more charged products, leading to more final 
fragments with a charge smaller than 5. Finally for all 
systems the detection reduces the fragment multiplicity 
by about 0.6 units. These final calculated multiplicities 
are in good agreement with the experimental values. The 
figure also shows that, although fixed at the dynamical 
stage, this observable cannot help in choosing the asys- 
tiffness of the EOS, because at each stage of the calcula- 
tion the fragment multiplicities are located on the same 
straight line whether the asysoft or the asystiff EOS was 
used. 



B. Preequilibrium emission 

We define the multifragmenting source as the largest 
single cluster recognized by the clusterisation algorithm 
at t=120 fm/c, the rest of the system being considered as 
preequilibrium emission. As in ref. [37] we find that the 
rate of nucleon emission in SMF is large between 60 and 
120 fm/c to become constant at a small value afterwards; 
this is a general feature of all semi-classical transport 
models. The evolution of the number of preequilibrium 
nucleons (defined as the difference between numbers of 
neutrons or protons of the system and of the source ) with 



FIG. 6. (color on line) Average neutron (circles) and pro- 
ton (stars) preequilibrium emission vs the N/Z of the to- 
tal systems at 45 ^4MeV. The two points corresponding to 
N/Z=1.38 have been shifted for visibility. Results are shown 
for the two asy-EOS. 



the N/Z of the systems is shown in Fig. 6, for reactions at 
45 A MeV. In all cases about 23% of the mass (charge) of 
the system has been emitted at 120 fm/c. The neutron- 
richest system has 24 extra-neutrons as compared to the 
neutron-poorest one; we observe that the former system 
ejected 8-9 more neutrons, and two protons less, With an 
asysoft EOS, there are more neutrons and less protons 
emitted than in the asystiff case. 

In the model most of that emission occurs during the 
expansion of the system, and thus it is a witness of sub- 
normal density EOS [37] ; in that density region the sym- 
metry potential is more repulsive for neutrons, and more 
attractive for protons in the asysoft case, which explains 
the relative values displayed in Fig. 6. Whatever the 
EOS is, in all cases the neutron to proton ratio of the 
fast emission is more neutron-rich than the system for 
136 Xe+ 124 Sn, and more neutron-poor for 124 Xe+ 112 Sn. 
Consequently the isospin content of the multifragment- 
ing sources is smaller than that of the total system for 
136 Xe+ 124 Sn, and larger for 124 Xe+ 112 Sn, which reduces 
the N/Z range explored in multifragmcntation. For the 
intermediate systems, the source is slightly more neutron- 
rich than the system with the asystiff EOS, and less for 
the asysoft EOS. 

Isolating preequilibrium emission in experiments is a 
difficult task. A possible choice is to look at high en- 
ergy nucleons, possibly including nucleons bound in clus- 
ters [6, 38]. We propose as an alternative to look at the 
complementary part of the proton preequilibrium emis- 
sion, the "liquid phase", through the value of the total 
charge bound in fragments, Zb 5 = \^* //ras Zz>5- Obvi- 
ously this variable is only meaningful provided that the 
detected events contain a large part of the total charge 
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TABLE III. For the 3 systems studied at E/A=45 MeV: mean experimental value of Zb 5 ; for SMF simulations, and the two 
asy-EOS: number of protons lost at 120 fm/c, Zb 5 for primary fragments (at 300 fm/c), for cold fragments and after filtering. 
In all cases the numbers between parentheses give the standard deviations of the corresponding distributions. 



System 


Exp. 


SMF 

asy (Z pr ) (Z p b :; m ) (Z' b ° ld } <Z£ W ) 


124 Xe+ 112 Sn 

136 Xe+ 112 Sn 
i36 Xe+ i24 Sn 


41.9 (8.1) 
44.4 (8.1) 
46.2 (8.3) 


stiff 24.60 (1.23) 63.2 (2.7) 47.2 (5.06) 41.6 (8.42) 
soft 23.72 (1.23) 63.9 (2.8) 47.4 (5.06) 41.6 (8.52) 
stiff 23.65 (1.21) 64.4 (2.8) 48.9 (5.21) 42.9 (8.72) 
soft 22.55 (1.12) 65.7 (2.6) 49.6 (5.17) 43.4 8.90) 
stiff 22.95 (1.20) 65.2 (2.9) 49.8 (5.48) 43.4 (8.96) 
soft 21.71 (1.13) 66.6 (2.7) 50.6 (5.42) 44.0 (9.06) 



of the system, which is the case of our selected events. 
We want to test the dependence of Zb 5 on the symmetry 
energy term implemented in the simulations. In this aim 
we have followed the evolution of Z^ as a function of the 
collision stage, for the measured systems, and compared 
the final result of the calculations with the experimental 
data. We show in Table III the mean values and standard 
deviations of the experimental and calculated (primary, 
and cold before and after filtering) distributions for the 
different systems. We clearly observe that the experimen- 
tal distributions arc shifted upwards for the neutron-rich 
systems. In simulations the mean values of the distribu- 
tions decrease whereas the standard deviations increase 
when going from primary to cold fragments and to the 
filtered distributions. The simulated mean values at all 
stages are, as in the experiment, larger when the system is 
more neutron-rich. The filtered mean value agrees rather 
well with the experimental one for 124 Xe+ 112 Sn, while it 
is more and more underestimated for the two other sys- 
tems. We also notice that the filtered distributions are 
broader than the measured ones. We remark that the fi- 
nal calculated (Z{, b ) does not depend on the asy-stiffness 
of the EOS for the neutron-poor system whereas it be- 
comes larger in the asy-soft case when the system con- 
tains more neutrons. 

If we compare the results for the extreme systems, 
i36 Xe+ i24 Sn and i24 Xc+ ii2 Sllj the difference between 

the number of prccquilibrium protons is negative (fourth 
column of table III and Fig. 6). Therefore we expect 
a positive difference, AZb 5 , when we turn to the values 
of Z bs , which is observed in the model for primary frag- 
ments at 300 fm/c, and for cold fragments. Note that the 
difference between the two systems grows larger during 
secondary decay, indicating that while it is essentially 
due to the first stage of the collision, part of it comes 
from the evaporation process. The filtering makes the 
AZb 5 smaller again. At all steps the difference between 
the two systems is larger for the asysoft EOS. Experi- 
mentally we do observe a positive AZ &5 , which is larger 
than what is obtained with any of the two EOS. The au- 
thors of ref [37] underline that the calculated AZb^ values 
depend on the fragment formation process; the spinodal 
decomposition which drives multifragmentation in SMF 
strongly disfavours the formation of light fragments as 



compared with an Antisymmctrizcd Molecular Dynam- 
ics calculation, AMD, (see figure 8 of ref. [37]). However 
in view of the significant measured difference between the 
two systems, we think that AZ b5 might be in the future 
a good observable for constraining the EOS. 



C. Fragment kinetic energies 

The calculated average kinetic energy of the final frag- 
ments at 45 A MeV is represented in Fig. 3. As in [32], 
we have added to the average kinetic energies a ther- 
mal term (3T/2) which is not contained in the calcula- 
tion. One immediately observes that the model largely 
underestimates the measured energies. The underestima- 
tion of the energy of fragments produced through multi- 
fragmentation in central collisions seems to be a general 
drawback of semi-classical transport models. In [32], we 
found that in the BOB simulation the fragment kinetic 
energies were about 20% below the experimental values 
for 129 Xe+" at Sn at 32 A MeV. Several causes of the dis- 
crepancy were identified [22]: thermal fluctuations were 
reduced by the large number of test particles, whereas 
quantal fluctuations were neglected. Second the finite 
range of the force introduced in the model was too large 
at normal density, creating surface energy at the detri- 
ment of fragment kinetic energy; the initial expansion 
energy was largely canceled. With the present SMF ver- 
sion the discrepancy between calculation and experiment 
reaches ^40%. In this model the range of the force is cor- 
rect for all densities; thermal fluctuations are still under- 
estimated and in addition the large amount of prccquilib- 
rium emission considerably reduces the available energy 
of the fragmenting system. Indeed in reference [37] it 
is shown that the AMD model produces less preequilib- 
rium nucleons, and fragments with larger kinetic ener- 
gies. The overestimation of particle emission seems to be 
a rather general problem of semi-classical approaches, as 
also stressed in [39]. However, it is interesting to note 
that in all simulations (BOB, SMF, AMD) for Sn+Sn or 
Xe+Sn collisions between 30 and 50 A MeV the excita- 
tion energy of the fragments is found close to 3 A MeV. 
The lack of available energy in SMF with respect to AMD 
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thus only appears in the kinetic energy of the fragments. 
Finally the isospin effect experimentally observed is not 
predicted by the calculation, which is not surprising if 
the effect comes from the expansion energy, which is con- 
sumed during the fragment formation. 

Coming back to the results displayed in table III, we 
notice that the values of Z b& from the filtered model and 
the experiment are similar, which would imply a cor- 
rect number of nucleons emitted at 300 fm/c. We should 
however stress again that we could not apply the exact 
experimental selections to the calculation (completeness 
of the events and use of a simple geometrical filter). 

V. SUMMARY AND CONCLUSIONS 

We have studied central collisions between different 
isotopes of Xe and Sn nuclei, leading to systems differ- 
ing by their number of neutrons, at incident energies of 
32 and 45 ^4MeV. Experimentally we highlighted sev- 
eral isospin effects. At a given incident energy the lcp 
multiplicity decreases whereas the fragment multiplicity 
increases with the N/Z of the system. The stochastic 
mean field model developed in Catania well accounts for 
the fragment partition properties at 45 AMeV, while at 
32 A MeV the model is not able to describe the multifrag- 
mentation process. Comparing model and experimental 
data, we infer that the increase of the fragment multi- 
plicity essentially comes from the dynamical part of the 
reaction, but does not help in constraining the symmetry 
energy term of the EOS. Indeed this seems to be related 
just to the larger size of the fragmenting source in the 
neutron-richer cases. 



We also found that on average the fragment kinetic 
energy grows larger when the isospin of the system in- 
creases, and we think that this observation confirms the 
existence of an expansion energy and shows that the pri- 
mary fragments keep some sign of the neutron-richness 
of the initial system. The model fails in reproducing the 
fragment kinetic properties, one of the reasons being that 
too much energy is removed from the system by preequi- 
librium emission. 

Finally we propose to use the charge bound in frag- 
ments as an alternate variable to get information on lcp 
preequilibrium emission. We experimentally observe a 
sizeable difference between the values of measured for 
136 Xe+ 124 Sn and 124 Xe+ 112 Sn, which makes it a promis- 
ing observable to constrain the symmetry energy. The 
calculated values of this difference depend on the asys- 
tiffness of the EOS, but are smaller than the measured 
one both in the asysoft and asystiff case. 

In prospective further developments of stochastic 
transport models are in progress. Phase space fluctu- 
ations are introduced by a stochastic treatment of the 
nucleon-nucleon collision integral. Once a two-body col- 
lision occurs, two clouds of test particles, which simu- 
late the wave-packet extension of the colliding nucleons, 
are moved simultaneously to new phase-space locations. 
This procedure, that corresponds to a numerical imple- 
mentation of the Boltzmann-Langevin equation, has been 
proven to give the correct amplitude of equilibrium ther- 
mal fluctuations [40]. First tests with this model show 
an improvement of the description of the kinematical and 
partition properties of fragmenting systems in dissipative 
collisions around the Fermi energy [41]. 
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